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^ The analysis of the Higgs boson data by the ATLAS and CMS Collaborations appears to exhibit 

an excess of h — > 77 events above the Standard Model (SM) expectations; whereas no significant 
excess is observed in h —> ZZ* — )> four lepton events, albeit with large statistical uncertainty due to 
the small data sample. These results (assuming they persist with further data) could be explained 
^s^J by a pair of nearly mass-degenerate scalars, one of which is a SM-like Higgs boson and the other 

is a scalar with suppressed couplings to and ZZ. In the two Higgs doublet model, the 

observed 77 and ZZ* — >> four lepton data can be reproduced by an approximately degenerate CP- 
O 1 even (h) and CP-odd (A) Higgs boson for values of sin(/3 — a) near unity and 0.7 <J tan/3 <J 1. An 

^ ^ enhanced 77 signal can also arise in cases where mu — mn, tuh — tua, or uiy h c± ttih — tua- Since 

qj the ZZ* — >• four lepton signal derives primarily from a SM-like Higgs boson whereas the 77 signal 

receives contributions from two (or more) nearly mass- degenerate states, one would expect a slightly 
different invariant mass peak in the ZZ* — > four lepton and 77 channels. The phenomenological 
£N| consequences of such models can be tested with additional Higgs data that will be collected at the 

LHC in the near future. 
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1. INTRODUCTION 



> 



The ATLAS pQ and CMS [2] collaborations have recently announced the discovery of a new boson with the properties 
approximating those of the Standard Model Higgs boson (h) [3]. The cleanest channels for the observation of the 
Higgs boson are via the rare decays h — » 77 and h — » ZZ* — >> £+£~£'+£'~ (where £ = e or /i). Evidence for Higgs 
boson production and decay is also seen in h — >> WW* — » £ + v£~V. Using this data, both ATLAS and CMS have 
compared the observed values of a x BR for the Higgs boson production and decay into the various final states, relative 
to the corresponding Standard Model predictions. Including both the 2011 and 2012 data sets, both collaboration 
observe excesses of the 77 channels [4j [5] , whereas the WW* and ZZ* channels are consistent with Standard Model 
expectations. 

Although the excess in the 77 channel is not yet statistically significant (at best it is a 2 sigma excess if the 
Higgs data from both experimental collaborations are combined), it is instructive to consider mechanisms that could 
enhance the 77 channel while not affecting the Standard Model-like results for pp into vector boson pairs. Note that 
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the Higgs boson couples to 77 through a one-loop process, whereas it couples to vector bosons at tree level. Thus, 
one could enhance the 77 signal by proposing new charged particles beyond the Standard Model that couple to the 
Higgs boson [6]. Such particles would contribute in the loop amplitude for h —> 77 and could potentially enhance 
the h — >• 77 branching ratio. However, in such a scenario one must check that the same particles do not contribute 
significantly to the gg — >• h production mechanism (which is also generated by a one-loop process). Otherwise, one 
might find that a x BR for Higgs production and decay into WW* and ZZ* is enhanced or suppressed with respect 
to the Standard Model in contradiction to the observed data. 

Alternatively, one could enhance the 77 signal by increasing the branching ratio for Higgs decay into 77 with 
respect to that of the Standard Model. One way to accomplish this is to reduce the corresponding Higgs branching 
ratios into other decay channels. The most effective mechanism of this type is one where the dominant Higgs decay 
rate into bb pairs is significantly reduced. Specific examples can be found in Ref. [7]. 

There exists a third mechanism that can potentially enhance the 77 signal while not affecting the WW* and ZZ* 
signal. For simplicity, we assume that the Higgs sector is CP-conserving. Suppose that there exist two nearly mass- 
degenerate scalars, a CP-even state h and a CP-odd state A. If the properties of h are approximately given by those 
of the Standard Model Higgs boson, then the observed a x BR(h —> VV*) [where V = W or Z] should be consistent 
with Standard Model Higgs production and decay, due to the fact that there are no tree- level AVV couplings. 1 On 
the other hand, the production cross section for gg —> A is typically larger than the production cross section for 
gg —> h, and the decay A — ^ 77 is mediated at one- loop by the top quark loop. Thus, if ~ m^, then both the h 
and the A will appear in the 77 signal, thereby yielding an enhanced 77 signal. 

It is also possible to achieve the same enhancement in a CP-conserving Higgs sector that only contains Higgs doublets 
(and possibly Higgs singlets) with no higher Higgs multiplets. If there exists a CP-even state h\ whose properties are 
approximately given by those of the Standard Model Higgs boson, then this one state will approximately saturate 
the sum rule for hVV couplings [9], in which case the couplings of any other neutral Higgs boson (both CP-even 
and CP-odd) to VV will be highly suppressed. In such a case, if two CP-even states h\ and h® are approximately 
mass-degenerate, then both can contribute to an enhanced 77 signal observed at the LHC, whereas only h\ production 
and decay yields a significant ZZ* — >• four leptons signal. 

In this paper, we consider the implications of the CP-conserving two Higgs doublet extension of the Standard 
Model (2HDM). There is already a substantial literature that analyzes the implications of the present Higgs data in 
the framework of the 2HDM [TOj [TTj . We shall address the question of whether the presence of nearly mass-degenerate 
neutral scalars of the 2HDM can be responsible for an enhanced 77 signal. Such a mechanism was first proposed in 
Ref. [12] in the context of the non-minimal supersymmetric extension of the Standard Model (NMSSM), and was also 
recently applied in the context of the 2HDM in an independent study [T3] . 

What is the origin of the approximate mass degeneracy of the two near-degenerate Higgs states? In most cases, the 
near-degeneracy is accidental. However, one could imagine that a near-mass-degenerate CP-even/CP-odd scalar pair 
might originate from a single complex scalar $>°. The presence of a global U(l) symmetry would then yield identical 
masses for the real and imaginary parts of 3>°. If the explicit breaking of the global U(l) symmetry is small (perhaps 
due to loop corrections), the resulting CP-even and CP-odd scalars would end up as approximate mass-degenerate 
states. Whether the near- degeneracy is accidental or the result of an approximate symmetry, it is important to analyze 
the LHC Higgs data to determine if such a scenario is realized in nature or can be excluded. In Ref. [14], a set of 
diagnostic tools was developed that can be used to determine whether mass-degenerate Higgs states are present in 
the LHC Higgs data. 

We begin our analysis by computing the a x BR for the production and decay of h and A separately, under the 
assumption that ~ tua — 125 GeV and the couplings of h to W + W~ and ZZ approximate the corresponding 
Standard Model Higgs boson couplings. We then determine the region of 2HDM parameter space in which the sum 



An interesting caveat, investigated in Ref. [8], is that there is a very large enhancement of A — y Z7, feeding into the A — y 4£ used 
experimentally to detect ZZ* . In that case, only A —y WW* would be suppressed. 
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of the h and A signals into 77 match the enhanced rates suggested by the central values measured by the ATLAS 
and CMS Collaborations. We also check to see whether an excess in the 77 channel can be achieved with other mass 
degenerate scalars (ft, H), (H, A) and (ft, H, A). Under the assumption that mass-degenerate neutral Higgs states 
can be invoked to explain the excess in the 77 channel, we predict which additional Higgs channels must also exhibit 
deviations from their corresponding Standard Model predictions. 



2. SETTING UP THE 2HDM PARAMETER SCAN 

The 2HDM consists of two hypercharge-one scalar doublet fields, denoted by <£ a , where a = 1, 2 [15 . The Lagrangian 
of the most general 2HDM contains a scalar potential, with two real and one complex scalar squared-mass term and 
four real and three complex dimensionless self-couplings, and the most general set of dimension-4 Higgs- fermion 
Yukawa interactions. For example, the most general Yukawa Lagrangian, in terms of the quark mass-eigenstate fields, 



J2V = Ul^Ur + D L K^-rf a U R + U L K<S>+ V a*D R + D L <S> aV ^D R + h.c. . 



(2.1) 



where <fr a = (3>° , <fr~) = icr 2 &* anc ^ K ^ s ^ ne CKM mixing matrix. In eq. (2.1), there is an implicit sum over the 
index a = 1, 2, and the rj u,D are 3x3 Yukawa coupling matrices. However, such models generically possess tree- level 
Higgs- mediated flavor- changing neutral currents (FCNCs), in conflict with experimental data that requires FCNC 
interactions to be significantly suppressed. 

To avoid dangerous FCNC interactions in a natural way, we impose a Z2 symmetry on the dimension-4 terms of the 



Higgs Lagrangian in such a way that removes two of the four Yukawa coupling matrices of eq. (2.1 ). The corresponding 



symmetry consists of choosing one of the two Higgs fields to be odd under the Z2 symmetry. The resulting Higgs 
potential is given by: 



V = m?!*!*! + m 2 22 $i<f>2 ~ (m? 2 *i*2 + h.c.) + |A X (V^i)' + \\ 2 ($^ 2 ) 2 



(2.2) 



where we have allowed for a softly-broken Z2 symmetry due to the presence of the term proportional to m\ 2 . For 
simplicity, we shall assume that the potentially complex terms in eq. (2.2), m\ 2 and A5 are real, in which case the 



Higgs potential is CP-conserving. In addition, some of the fermion fields may also be odd under the Z2 symmetry, 
depending on which two of the four Yukawa coupling matrices are required to vanish. Different choices lead to different 
Yukawa interactions. In this paper, we focus on two different model choices, known in the literature as Type-I [T6l[T7] 
and Type-II [13 [18]. In the Type-I 2HDM, tjY 
exclusively to $ 2 . In the Type-II 2HDM, rjf 



in eq. (2.1), which means that all quarks and leptons couple 



in eq. (2.1), which means that the up-type quarks couple 



exclusively to $2 and the down- type quarks and charged leptons couple exclusively to <3>i. 

We assume that the parameters of the Higgs potential are chosen such that the SU(2)xU(l) electroweak symmetry 
is broken to U(l) electromagnetism. The neutral Higgs fields then acquire vacuum expectation values, = V a/V2, 
where v 2 +v 2 — 4m^ / g 2 = (246 GeV) 2 and tan/3 = v 2 jv\. By convention, we can choose < /3 < \k. Expanding the 
Higgs potential about its minimum, we then diagonalize the resulting scalar squared-mass matrices. Three Goldstone 
boson states are eaten by the and Z gauge bosons, leaving five physical degrees of freedom: a charged Higgs 
pair, i^, two CP-even neutral Higgs states, h and H (defined such that < ra#), and one CP-odd neutral Higgs 
boson A. In diagonalizing the CP-even neutral Higgs squared-mass matrix, one also obtains a CP-even Higgs mixing 
angle, a. By convention, we take \a\ <tt/2. In the Type-I and Type-II 2HDM, the Higgs- fermion couplings are flavor 
diagonal and depend on the two angles a and (3 as shown in Table [TJ 

The aim of this paper is to scan the Type-I and Type-II 2HDM parameter spaces allowing the 77 signal resulting 
from the production of two near-mass-degenerate Higgs states to be enhanced above the Standard Model (SM) 
rate, while assuming that the corresponding ZZ* — >• 4 leptons (and WW*) signal is approximately given by the 
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Type-I 
h A H 


Type-II 
ft AH 


Up-type quarks 
Down-type quarks and charged leptons 


cos a 1 sin f3 cot /3 sin a / sin /3 
cos a 1 sin /3 — cot /3 sin a / sin /3 


cos a 1 sin /3 cot f3 sin a / sin (3 
— sin a/ cos /3 tan/3 cos a/ cos f5 



TABLE I: Couplings of the fermions to the lighter and heavier CP-even scalars (ft and H), and the CP-odd scalar (A). 



corresponding SM rate. In general, given a final state /, we define Rf as the ratio of the number of events predicted 
in the 2HDM with near-mass-degenerate Higgs states to that obtained in the SM: 



Rf 



N(pp /) 2 HDM 



N(pp /) SM 

Denoting the two near-mass-degenerate Higgs states by hi and ft,2, we obtain 



Rf = R 1 } 1 +R h f 2 



where 



R 



f 



a(pp -)> S) 2 udm BR(S -)> /) 2 HDM 
a(pp -> /i S m) BR(ft S M -> f) 



(2.3) 



(2.4) 



(2.5) 



for S = hi , h2 and a is the Higgs production cross section, BR the branching ratio, and hsM is the SM Higgs boson. 
In our analysis, we include all Higgs production mechanisms, namely, gluon-gluon fusion using HIGLU at NLO [19], 
vector boson fusion (VBF) [20], Higgs production in association with either W, Z or ti [20], and bb fusion [21]. 

Since we do not expect Higgs coupling measurements at the LHC to be better than about 20% after the full 2012 
data set is analyzed, we have performed our 2HDM scans under the assumption that 



0.8 < R ZZ < 1.2. 



(2.6) 



In the case where hi = h and h 2 = A the constraint imposed by our Rzz assumption implies that the ZZ* — > 4 leptons 
signal is due almost entirely to the production and decay of ft, due to the absence of a tree-level AZZ coupling. That 
is, the properties of the ft are SM-like. Remarkably, the same conclusions emerge in the case where both hi and hi are 
CP-even states. In this case, we find that one of the two CP-even states is SM-like, where the other has suppressed 
tree-level couplings to ZZ. 

We perform our 2HDM scans separately for the Type-I and Type-II Higgs- fermion couplings, by scanning over the 
neutral and charged Higgs masses, the angles a and j3 and the soft-breaking parameter m^ 2 , subject to the constraint 



of eq. (2.6) and the near mass-degeneracy of two (or three) neutral Higgs bosons. If the two nearly mass-degenerate 
Higgs states are CP-even, then we shall assume that the mass difference of these states is large enough (compared 
to the corresponding Higgs boson widths) so that we can neglect possible interference effects in the production and 
decay process. 2 This is not a significant constraint on our analysis, as the Higgs widths are significantly smaller than 
the experimental mass resolution of the ATLAS and CMS experiments. 

In addition, the allowed points of our scans must also satisfy the known indirect experimental constraints on the 
Type-I and Type-II 2HDM. These constraints are described in Section [3] After imposing these requirements, we 
identify the surviving parameter regimes where the 77 signal is enhanced, and discuss additional Higgs phenomena 
that must be observed in order to confirm the scenarios advocated in this paper. 



In the case of a mass-degenerate CP-even/CP-odd Higgs pair, there is no interference due to CP invariance. 
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3. CONSTRAINTS ON THE 2HDM PARAMETERS 

In this section, we summarize the indirect experimental constraints on the most general CP-conserving 2HDM's 
with Type-I and Type-II Higgs-fermion Yukawa couplings, respectively. First, we scan over the parameter space 
subject to the following three constraints: the Higgs scalar potential is (i) bounded from below [22]; (ii) satisfies 
tree- level unitarity [23 j; and (iii) is consistent with constraints from the Peskin-Takeuchi S and T parameters [24^ [25] 
as derived from electroweak precision observables [26-28 ]. Additional indirect 2HDM constraints arise from charged 
Higgs exchange contributions to processes involving the b quark and r-lepton. Such processes yield constraints in 
the m#±-tan/3 plane. Consequently, if the rates for these processes are consistent with SM expectations, then the 
corresponding 2HDM constraints are (typically) relaxed in the limit of large m H ± due to the decoupling properties of 
the 2HDM [29]. At present, there is an anomaly observed by the BaBar collaboration in the rate for B — >> D^r~V T 
which deviates by 3.4 a (when D and D* final states are combined) form the SM prediction [30]. The observed 
deviation cannot be explained by a charged Higgs boson in the Type-II 2HDM. Indeed the analysis of Ref. [30] 
excludes the Type-II 2HDM for any value of tan/3/m#± at the 99.8% CL. Until this observation is independently 
confirmed by the BELLE collaboration, we will not include this result among our 2HDM constraints. 

A number of other observables in B physics provide constraints in the m#±-tan/3 parameter plane. The most 
precise SM prediction for = T(Z — » bb)/T(Z — >• hadrons), including electroweak two-loop and QCD three-loop 
corrections, deviates by two standard deviations from the experimentally measured value [31 . The inclusion of 
2HDM contributions could alleviate this discrepancy in particular regions of the parameter space [32 -35 . Another 
constraining observable is the decay rate for b — » sj. The theoretical prediction for b — >• sj in the SM has been 
performed to order 0(a 2 s ) [36] and is in good agreement with the experimental result [37]. However, due to theoretical 
uncertainties coming mainly from higher order QCD corrections and the uncertainties in the CKM matrix element, 
there is still some room for 2HDM contributions. Other B observables such as B-B mixing data and the rate for 
B + — >• t + v t also provide constraints on the 2HDM parameter space. 

In the Type-II 2HDM, provides the most restrictive bound in the ra#±-tan/3 plane in the small tan/3 region, 
as shown in Refs. [27] [35]. However, the strongest bound on the charged Higgs mass derives from the b — > 
measurement, which yields m H ± > 340 GeV almost independently of the value of tan/3 [38]. In contrast, B-B mixing 
yields less constraining bounds in the m#±-tan/3 plane, whereas B + — >• t + v t constrains mainly large values of tan/3 
for small charged Higgs masses [39 . In the Type-I 2HDM, the most restrictive bound in the ra#±-tan/3 plane is due 
to b 87 [39 . In contrast to the Type-II model, there is no strong bound (independent of the value of tan/3) on 
the charged Higgs mass. Finally we note that all other experimental bounds on the Type-I 2HDM from 5-physics 
observables were shown to be less restrictive [27] [39] . 

With the exception of B — >> t + v t and B — > D^r~V T (where charged Higgs exchange contributes at tree- level), 
the 2HDM contributions to B physics observables arise via one-loop radiative corrections. For such observables, it 
is always possible that constraints on the 2HDM parameter space could be relaxed due to cancellations in the loop 
from other sources of new physics. One of the well-known examples of this phenomenon is the partial cancellation 
of the charged Higgs loop and the chargino loop contributions to b — » in the MSSM [40]. Thus, to be flexible in 
our presentation, we shall present results of our 2HDM scans with and without the bounds from 5-physics and R^. 
In the former case, the 95% CL bounds from 5-physics and R^ are employed. We will denote our parameter scan 
as "constrained," when all experimental and theoretical bounds are considered, and "unconstrained" when all except 
the Rb and -physics bounds are taken into account. 



DEGENERATE h AND A 



We begin by fixing tua = rrth = 125 GeV. We scan over 2HDM parameters subject to eq. (2.6). As a consequence, 
the h couplings to ZZ and ti cannot differ much from their SM values. The 77 signal is enhanced if the 

contribution from the A is large; this can be achieved for values of tan/3 < 1. Thus, we focus our analysis on 
the region 0.5 < tan/3 < 2. Since A does not couple to ZZ, the square of the hZZ coupling, i.e. sin 2 (/3 — a), is 
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constrained to be near its SM value. In the tan/3 parameter regime of interest, we find that 0.7 < sin 2 (/? — a) < 1 
for both Type-I and Type-II Higgs-fermion Yukawa couplings. 3 

By convention (and without loss of generality), we define a and f3 such that — \i\ < a < \i\ and < f3 < Thus, 
for 0.5 < tan/3 < 2 it follows that 0.84 < sin(/3 — a) < 1. In particular, negative values of sin(/3 — a) near —1 are 
not permitted. In this parameter regime, | cos(/3 — a)\ < 0.55 with both signs allowed. Given these constraints, the 
range of possible values of a is also constrained. Using the trigonometric identities, 

sin a = sin /3 cos(/3 — a) — cos /3 sin(/3 — a) , (4.1) 
cos a = cos j3 cos(/3 — a) + sin f3 sin(/3 — a) , (4.2) 

one can easily determine the allowed range of possible a values in our scans. 

In our 2HDM parameter scan, the mass of the heavier CP-even scalar (ra#) is kept between 200 and 1000 GeV, and 
a is allowed to vary subject to the constraints implicit in eqs. (4.1) and ( |4.2| ). As discussed in Section [3j we impose 



the requirements of a scalar potential that is bounded from below and satisfies unitarity, and apply the constraints 
from precision electroweak observables. We also vary the charged Higgs mass between 500 and 1000 GeV subject 
to all the constraints discussed in Section [3] After imposing all the relevant constraints, we note that the one-loop 
diagrams mediated by the charged Higgs boson contribute very little to pp — >• (h,A) —^77. 



4.1. Type-I 2HDM 

In the simulation of a mass-degenerate h, A pair in the Type-I 2HDM, we find sin (f3 — a) > 0.88 for tan (3 ~ 0.8 and 
sin (j3 — a) > 0.93 for tan/3 ~ 1.6. These results are easy to understand in the approximation where Higgs production 
is due exclusively to gluon-gluon fusion via the top quark loop, and the total Higgs width is well approximated by 
r(ft — » 66). In this simplified scenario, 

„ cos 2 a . 9 , _ » sin 2 B . 9 , ^ N , ^ 

Rzz = —Y-r sin 2 (a - /?) = sin 2 (a - P . 4.3 

sur f3 cos 2 a 

Using the Type-I couplings of Table |IJ the first factor derives from the square of the tih coupling that governs the 
gluon-gluon fusion cross section, the second from T(ft — ^ ZZ) and the third from r to tai ~ r(ft —> 66) which appears in 



the denominator. Thus sin 2 (a — f3) = Rzz > 0.8 as a result of the constraint imposed by eq. (2.6). This precludes 
small values for sin 2 (a — and we expect SM-like couplings of ft to gauge bosons. 

We have generated over 7000 points in parameter space satisfying all constraints. The results for R 11 as a function 
of tan/3 are presented in Fig. [l] for the unconstrained scenario (left panel) and the constrained scenario (right panel). 
We note that R^ lies below unity regardless of tan /3. This result was previously obtained in Ref. [10] in the limit 
of vfiA ^> Triht and in Ref. [41], allowing for the mixing with the heavier CP-even and CP-odd scalars. In contrast, 
R^ increases as tan/3 decreases, as expected, due to the increase in the Att coupling [cf. Table [ij. As a result, the 
observed rate for pp — ^ (ft, A) — » 77 normalized to the corresponding SM rate can take values as large as 2, for 
tan /3 ~ 0.6. However, constraints from 5-physics restrict the normalized rate for pp — >• (ft, A) — » 77 to a maximal 
value of approximately 1.4, for tan/3 ~ 1. 

Since A does not couple to VV, the 77 decays detected in VBF production can only be due to the ft intermediate 
state. The ATLAS and CMS experiments can isolate Higgs signal events with a set of additional criteria (e.g. events 
with two forward jets with certain transverse momentum cuts and a central jet veto) which they designate as VBF 
Higgs events. In practice the experimental VBF Higgs events have significant contamination 4 of Higgs events produced 



3 This r esult is not surprising. For values of tan (3 ~ 1, the production cross-section for h production is SM-like. In light of the constraint 
of eq. (|2.6fc, it follows that the hZZ coupling should be approximately SM-like, which implies a value of sin(/3 — a) close to 1 [3 29 . 

4 The ATLAS and CMS collaborations quote contamination rates for the gluon-gluon fusion Higgs events of roughly 30% ^ [2] , although 
in practice this number has a rather large error bar. 
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FIG. 1: Left panel: i? 77 as a function of tan/3 for h (blue/black), A (green/gray), and the total observable rate (cyan/light- 
gray), obtained by summing the rates with intermediate h and A, for the unconstrained scenario. Right panel: Total rate for 
R 11 as a function of tan/3 for the constrained (red/black) and unconstrained (green/gray) scenarios. 



by gluon-gluon fusion with the subsequent radiation of two additional jets. Nevertheless, in this work we find it 
convenient to define a theoretical quantity, 

dVBF _ a (PP ~> VV ^)2HDM BR(fe -)> 77) 2 HDM ( . 

77 ~ <r(pp ^VV^ h SM ) BR(/i SM -+ 77) ' v 



which would be appropriate if VBF Higgs events could be identified with no contamination. This will prove sufficient 
for our purposes in this initial study. Likewise, we shall denote R 71 = R^ + R^ following the definitions given in 
eq. ([275b. 

Fig. [2] shows the allowed region in the i?^ F -i? 77 plane, for our set of points, for the constrained (red/black) and 
unconstrained (green/gray) scenarios. In the unconstrained scenario, contrary to what happens for the total i? 77 , the 
h intermediate state can induce values of R^ F larger than unity. The dispersion of points shows that, in the region 



Model I, h and A degenerate. Green - unconstrained; red - constrained 
2.2 I 1 1 1 




1 

D VBF 



FIG. 2: Allowed region in the i? 77 3F — R 11 plane with (red/black) and without (green/gray) the B-physics constraints. 
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of parameter space that we have studied, a comparison between R^ F and i? 77 is unlikely to allow an exclusion of the 
Type-I 2HDM. Indeed, for R 71 ~ 1.5, this model allows for any value of R^ F between 0.2 and 1.7. Only if R 71 ~ 2.0 
can we exclude values of R^ F above unity. 

In the constrained scenario, R^ F is now between 0.7 and 1.6, while R 77 lies in the range from 1 to 1.4. The allowed 
values (red/black) clearly show that an enhancement in i? 77 can only be achieved for R^ F close to 1. Conversely, 
large values of R^ F are attained only for a SM-like R 11 . Naively, it seems puzzling that values of i^X? F above 1 
are possible. After all, this quantity is only sensitive to h production and decay, and we are assuming that the hVV 
coupling is close to its SM value. A closer examination of the parameter scan reveals that the range of allowed a is 
rather limited, —0.9 < sin a < —0.4. For the values of sin a close to —1, the Type-I couplings of h to fermion pairs 
is suppressed relative to the corresponding SM couplings. Consequently, the partial width of h — > bb is reduced and 
so the corresponding branching ratio for h — ^ 77 is enhanced relative to its SM value. 5 



Model I, h and A degenerate. Green - unconstrained; red - constrained Model I, h and A degenerate. Green - unconstrained; red - constrained 



DC 




0.8 1 1.2 1.4 1.6 1.8 2 2.2 0.8 1 1.2 1.4 1.6 1.8 2 2.2 



FIG. 3: Left panel: Total R TT (h and A summed) as a function of i? 77 for the constrained (red/black) and unconstrained 
(green/gray) scenarios. QCD corrections aside, this figure also holds approximately for Rbb- Right panel: R^ b u (h and A 
summed) as a function of R 11 for the constrained (red/black) and unconstrained (green/gray) scenarios. 

The left panel of Fig. [3] exhibits our results for the inclusive rr final state, summing over all production mechanisms, 
for the constrained (red/black) and unconstrained (green/gray) scenarios. Note that both R TT and R 11 include 
contributions from h and A. For the unconstrained scenario, we see that if R 11 ~ 1.5 then the total R TT contribution 
cannot be smaller than about 3.5 and can be as large as 8. When the -physics constraints are included the main 
difference is again that R 11 < 1.4, whereas R TT lies in a very narrow band heavily dependent on the particular value 
of i? 77 but always below about 3. The correlation between the enhanced 77 and rr signals is a noteworthy prediction 
for this scenario. 

There are only very slight differences between R TT and R^, due to QCD corrections that distinguish the two 
processes, so the left panel of Fig. [3] would apply as well to a hypothetical measurement of the inclusive bb final state. 
Unfortunately, the detection of h — >• bb via gluon-gluon production at the LHC is swamped by background and is not 
possible in the inclusive mode. However, both the ATLAS and CMS collaborations expect some sensitivity to the bb 
final state in the production of the Higgs boson in association with a or Z (where leptonic decays of the vector 
bosons can be used to tag the event). The right panel of Fig. pi exhibits our results for the bb final state obtained 



5 There is a small additional enhancement to the /177 partial width since the contribution of the top quark loop (which negatively interferes 
with the dominant loop) is also reduced. 



8 



4.1 Type-1 2HDM 



4 DEGENERATE h AND A 



through VH production, ii^ H , with (red/black) and without (green/gray) the 5-physics constraints. Since A is not 
produced by this mechanism, we only get a contribution from the SM-like h and the enhancement in the bb channel 
disappears. Furthermore, in the constrained scenario, even an enhancement in R 71 corresponds to below 1. 



4.2. Type-II 2HDM 

In the analysis of a mass-degenerate h, A pair in the Type-II 2HDM, we have generated close to 1500 points. We 



again anticipate that sin(/3 — a) should be near 1. Applying the same simplified scenario that yielded eq. (4.3) for the 
Type-I 2HDM, we now obtain 



Rzz 



cos 2 a 
sin 2 /3 



sin (a — /3) — 



cos 2 /3 



sin 2 a 



(4.5) 



where we have employed the Type-II couplings of Table [TJ Assuming that cos(/3 — a) is small, we can expand in this 
small quantity by making use of 



cos 2 a cos 2 /3 
sin 2 /3 sin 2 a 



sin(/3 — a) + cot /3 cos(/3 — a) 
sin(/3 — a) — tan /3 cos(/3 — a) 



4 cos(/3 — a) 
sin 2/3 



+ (9(cos 2 (/3-a)) 



Inserting this result into eq. (4.5) yields 



cos(/3-a) ~ Um2f3(R zz -l). 



(4.6) 



(4.7) 



In light of eq. (2.6), the assumption that cos(/3 — a) is small is justified. In our simulation we find that the constraint 



imposed by eq. (2.6) leads to sin (/? — a) > 0.996, which is an even tighter restriction than suggested by eq. (4.7). 

Fig. [4] shows i? 77 as a function of tan /3 for the Type-II 2HDM. Here the h contribution by itself can only reach 
unity and, as in the Type-I 2HDM, the contribution from A becomes dominant for low tan f3 values. The total R 11 
can be as large as 3.5 for tan/3 ~ 0.7. Moreover, the requirement that 0.8 < Rzz < 1-2 leads to the exclusion of any 
point with tan/3 < 0.7 in our simulation. Because the most important constraint for the Type-II 2HDM is the one 
from Rb and since the charged Higgs mass is varied from 500 to 1000 GeV, tan ft will not be affected by any further 
constraints form B physics. This means that all results presented for this Type-II scenario already incorporate all 
available experimental and theoretical constraints. 



Model II, h and A degenerate. Blue - h; green - A; cyan - h + A 





i i i i 
















^ ur . v «„.^« ... . 


i i i i i i i 1 ■ - 



1.1 

tanp 



FIG. 4: R 11 as a function of tan/3 for h (blue/black), A (green/gray), and the total observable rate, obtained by summing the 
rates with intermediate h and A (cyan/light-gray). 
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Model II, h and A degenerate 



i 1 1 1 1 r~ 




0.82 0.84 0.86 0.88 0.9 0.92 0.94 0.96 0.£ 
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4 
3.5 
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2.5 
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Model II, h and A degenerate. Blue - h; green - A; cyan - h + A 




2 
R 



2.5 



3.5 



FIG. 5: Left panel: Allowed region in the R^^ F — R 11 . Right panel: R TT as a function of R 11 for ft (blue/black), A (green/gray), 
and the total observable rate, obtained by summing the rates with intermediate ft and A (cyan/light-gray). 



The left panel of Fig. [5] exhibits the allowed region in the R^ F -R 77 plane, for our set of points. In stark contrast 
with the Type-I 2HDM, here R^ F is predicted to lie in a rather narrow region close to 0.9. This is due to the much 
narrower range obtained for sin (/? — a), which is constrained to be near 1 in the Type-II scenario, and implies that 
cos a ~ sin f3 and sin a ~ — cos f3 (i.e. the decoupling limit). Consequently, the Type-II couplings of ft to fermion pairs 
are close to their SM values. It immediately follows that i?X? F should be close to 1 as indicated by the left panel of 
in Fig. [5] If the experimental data were to indicate a similar and sizeable enhancement of both R^ F and i? 77 , then 
the nth — tua Type-II scenario would be excluded. 

As for i? rr , the largest allowed values are less extreme than those of the Type-I model, as can be seen from the 
right panel of Fig. [5] The lower limit on R TT depends on the value of i? 77 . For example, for i? 77 ~ 1.5, the total R TT 
cannot be smaller than about 3.5. Notice that the right panel of Fig. [5] shows a rather strong correlation between 
R TT and i? 77 in the Type-II model, while the Type-I model points exhibited in the left panel of Fig. [3] are much more 
dispersed in the R 11 -R TT plane. 



5. DEGENERATE ft AND H 

We now turn to the possibility that rrih — run- We have generated sets of parameters such that tua and ra#± lie 
above 500 GeV. We have found that the 5, T, and U constraints force m H ± ~ m^, to within 10%. In this case, an 
enhancement of the 77 signal rate is visible (although not as pronounced) even for values of tan j3 somewhat larger 
than 2. As a result, we have focused our scan in a tan/3 regime between 0.5 and 5. 



5.1. Type-I 2HDM 

We begin with the Type-I 2HDM. In the case of a mass degenerate ft, A pair discussed previously sin (/? — a) > 0.85 
in both the Type-I and Type-II scenarios, with only a mild correlation with tan/3. In contrast, employing the same 



simplified scenario that was used in obtaining eq. (4.3), we now obtain 



^ h ^ H cos 2 a . 2 /n N sin 2 /3 sin 2 a 2 . sin 2 f3 . . . 

RZZ = Rzz + RzZ = —27i Sm (P - O)^^ + COS 2 (/? - a) = 1 ? ( 5 .1 ) 

sin p cos z a sin p sin a 
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Model I, h and H degenerate. Green - unconstrained; red - constrained Model I, h and H degenerate. Green - unconstrained; red - constrained 




FIG. 6: Left panel: Allowed region in the i? 77 3F — i? 77 plane for the constrained (red/black) and unconstrained (green/gray) 
scenarios. Right panel: Allowed region in the R 11 — R TT plane for the constrained (red/black) and unconstrained (green/gray) 
scenarios. 



where we have used the couplings of Table [T[ In fact, even using all production mechanisms, we find that Rzz ~ 1- 
Consequently, eq. ( |2.6| ) places almost no constraint on the value of sin (/? — a). 

The left panel of Fig. [ej exhibits the allowed region in the R^ F -R 71 plane, for the constrained (red/black) and 
unconstrained (green/gray) scenarios. The total R 11 cannot differ significantly from its SM value. In particular, the 
maximal allowed 77 enhancement is about 1.3 for the unconstrained scenario, and none of the enhanced values survive 
after imposing the constraints from B physics. In contrast, R^ F can be either very close to vanishing or take values 
as large as 4 in the unconstrained scenario and as large as 3 in the constrained scenario. 

The possibility of an enhanced R^ F arises due to the fact that a is not especially constrained by the requirement 
of eq. (2.6). In particular, the limit of sin a = ±1 [sin a = 0] corresponds to the fermiophobic limit for h [H]. In 
this limit, the dominant fermiophobic Higgs boson decay channels are WW* and ZZ*, in which case we expect the 
corresponding branching ratio into 77 (which should be approximately given by the ratio of the 77 and WW* partial 
widths) to be enhanced by a factor of roughly 5 relative to its value in the SM. Thus, even though the VBF production 
cross-section is roughly given by its SM value, it is not surprising that one can achieve values of R^ F as large as 
shown in the left panel of Fig. [6] Moreover, one also expects a reduced value of i? 77 , since in the fermiophobic 
limit the size of the gluon-gluon fusion cross section (which depends on the coupling of the Higgs boson to tt) is 
significantly suppressed. Assuming that the VBF production mechanism now dominates, it follows that the Higgs 
production cross-section has been reduced by roughly a factor of 10 relative to the SM Higgs production cross-section 
via gluon-gluon fusion. This reduction factor of 10 cannot be completely compensated by an increase in the Higgs to 
77 branching ratio, which implies that i? 77 < 1 in the fermiophobic regime, as shown in the left panel of Fig. [6] 

The right panel of Fig. [6] exhibits R TT as a function of R 11 . It is interesting to compare the Type-I scenarios 
exhibited by the left panel of Fig. [3j which holds for the case of a mass-degenerate h, A pair, with the right panel 
of Fig. |6j which holds for the case of a mass-degenerate h, H pair. In contrast to the former scenario, where i? 77 
could reach about 2 while R TT could reach 11, in the present scenario R 11 can be at most 1.3 while R TT is smaller 
than 5.5. Nevertheless, the lower bound on R TT for i? 77 > 1.1 is more stringent in the latter scenario than in the 
former. Indeed, when rrih — tua and R 11 = 1.2, we predict R TT > 2, whereas when rri} h ~ run and R 11 = 1.2, we 
predict R TT > 3.5. Again, in the constrained scenario, R 77 has to be below 1 while R TT cannot be above 2. The 
constrained scenario makes a very strong prediction: i? 77 has to be SM-like or smaller, whereas R T 
the SM prediction. 
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5.2. Type-II 2HDM 



In this case there is a strong correlation between sin (/3 — a) and tan j3. Employing the same simplified scenario 
that was used in obtaining eq. (4.3) and using the couplings of Table [IJ we now obtain 



R 



zz 



n ZZ 



Rzz 



■ sin 2 (/? — a) 



cos 2 (j 
sin 2 a 



sin 2 a o / n \ cos2 P 
■ cos (p — a) = — 



• 2 q °V . 2 1 • 2 q 

sin p sm a sin p 

The above result can be expressed directly in terms of tan/3 and sin(/3 — a) by employing eqs. (4.1) and (4.2), 

[tan 2 f3 + cot 2 /3 



2] C ; 



[t&n f3(cj_ a - sj_ 
sin(/3 



+ [tan/?(c|_ a - s|_J - 2c P - a s p . 
.J - (1 - tan 2 j3)c!3- a Sf}- a } 2 



(5.2) 



(5.3) 



where cp- a = cos(/3 — a) and = sin(/3 — a). Thus, eqs. (2.6) and ( |5.3[ ) impose a complicated constraint in the 

sin(/3 — a)-tan/3 plane. For example, for model points where 0.8 < R 11 < 1.5, the allowed regions in the sin(/3 — a)- 
tan /3 plane are shown in the left panel of Fig. [7j The uppermost (lowermost) two bands correspond to a dominant h 
(H) contribution to i? 77 . In the band around sin (/3 — a) ~ 1, h is gaugephilic and is gaugephobic. In the second 
band below, the contribution of h to i? 77 is still dominant, with a small contribution by H. In the next band, the 
roles of h and H are reversed, and in the band around sin (/3 — a) ~ 0, H is gaugephilic while h is gaugephobic. 




FIG. 7: Left panel: Values obtained in the tan/3 — sin (f3 — a) plane for the points generated, which satisfy 0.8 < R 11 < 1.5. 
Right panel: Values for R 11 as a function of tan/3 for the constrained (red/black) and unconstrained (green/gray) scenarios. 

The band structure exhibited in the left panel of Fig. [7] impacts the R 77 dependence on tan /3, as shown in the right 
panel of Fig. [TJ The generated set of points are plotted in the R^ F -R 77 plane in the left panel of Fig. [8j whereas 
the allowed regions in the i? 77 - J R rr plane is shown in the right panel of Fig. [8] for the constrained (red/black) and 
unconstrained (green/gray) scenarios. Focusing on the Type-II 2HDM, it is instructive to compare the scenarios with 
rrih ~ rriA and with rrih — run, respectively. Comparing left panels of Fig. [5] and Fig. [8j we see that the former 
scenario allows for larger values of i? 77 . In contrast, the values of R^ F are very constrained in the former scenario, 
while in the latter scenario these values can range from to 1.1. Comparing the right panels of Fig. [5] and Fig. |§J we 
see that 2.5 < R TT < 4.5 in the case of rrih = mA, whereas 0.2 < R TT < 2.7 in the case of = m#. 

Once the constraints of B physics are applied, the most striking difference from the results presented above is that 
the value of R 11 in the constrained scenario must be below 1.5. The allowed ranges of R^ F and R TT are quite 
narrow. However, a definite measurement of R 71 = 1.5 would force R TT « 2.5. It is also interesting to compare 
the rrih — rriH Type-II model analyzed here with the rrih — niH Type-I model examined in the previous subsection. 
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Model II, h and H degenerate. Green - unconstrained; red - constrained Model II, h and H degenerate. Green - unconstrained; red - constrained 
3 1 1 1 1 1 1 1 1 1 3 I 1 1 1 1 1 




FIG. 8: Left panel: Allowed region in the i? 77 3F -i? 77 plane for the constrained (red/black) and unconstrained (green/gray) 
scenarios. Right panel: Allowed region in the R 11 -R TT plane for the constrained (red/black) and unconstrained (green/gray) 
scenarios. 

Here the enhancement in R 11 can be larger, whereas the upper bounds on the enhancements in R^ F and R TT are 
somewhat reduced. In summary, it is more difficult to generate an enhancement of i? 77 in the rrih ~ m# case as 
compared to the rrih ~ vtia case. On the other hand, the corresponding enhancement of R rr is also reduced. 



6. DEGENERATE H AND A 

Now we turn to the case where rrih < rnn — tua — 125 GeV. In the previous cases we have studied, we could 
keep m H ± > 500 GeV. In this case, there are no points with m H ± > 200 GeV that survive the constraints coming 
from boundedness from below, unitarity, and precision electroweak measurements. Since the constraints from b — >• 
imply m H ± > 340 GeV in the Type-II 2HDM, the whole of parameter space with rrin ~ vtia — 125 GeV is excluded 
in the Type-II model. Thus we examine the rriA — ttih scenario in Type-I 2HDM, keeping 0.5 < tan /3 < 5. 

If rrih < ra#/2, then the branching ratio for H — ^ hh would be close to one, greatly suppressing the ZZ* signal 
(to which A cannot contribute). As a result, we shall assume that rrih > run/^- In particular, we examine the 
mass range 65 GeV < rrih < 110 GeV. Such a low-mass h could have been produced at LEP via the ZZh vertex, 
which is suppressed by a factor of sin (/? — a) in the 2HDM. The model points surviving the LEP constraints [42] , 
in which eq. ( |2.6| ) is satisfied due to H production and decay, are shown in the tan/3-sin (/? — a) plane of Fig. [9] for 
the constrained (red/black) and unconstrained (green/gray) scenarios. Note that because the charged Higgs boson in 
now forced to be light, the 5-physics constraints imply a value of tan/3 above 2. As expected, the points are centered 
around sin (/3 — a) = 0, but values as large as | sin (/? — a)\ ~ 0.5 are possible. 

The generated set of points are plotted in the i?XP F -.R 77 plane in the left panel of Fig. 
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for the constrained 

(red/black) and unconstrained (green/gray) scenarios. Contrary to the previous scenarios with ~ tua or rrih — mn 
in the Type-I model, both R^ F and i? 77 now have a much wider range of variation and can reach 2.3 simultaneously. 
However, in the constrained scenario the R^ F values drop to a maximum of around 1.2 while i? 77 drops to a maximum 



close to 1. The allowed region in the R^-R rr plane for this set of points is shown in the right panel of Fig. 10 We 
note that the constrained scenario lives in a region very close to the SM prediction. 
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FIG. 9: Values obtained in the tan /3-sin (f3 — a) plane for the points generated for the constrained (red/black) and unconstrained 
(green/gray) scenarios. 



Model I, H and A degenerate. Green - unconstrained; red - constrained Model I, H and A degenerate. Green - unconstrained; red - constrained 
3 1 1 1 1 1 1 1 5 i 1 1 1 1 1 




1.5 
R 



FIG. 10: Left panel: Allowed region in the i?^ F -i? 77 plane for the constrained (red/black) and unconstrained (green/gray) 
scenarios. Right panel: Allowed region in the R 11 -R TT plane for the constrained (red/black) and unconstrained (green/gray) 
scenarios. 



7. DEGENERATE h, H AND A 

In this section we discuss the case where ~ tuh — tua — 125 GeV, i.e. all three neutral scalars are nearly 
mass-degenerate. As in the previous scenario with ra# ~ m^, the electroweak precision constraints force the charged 
Higgs to be light (below approximately 200 GeV). Hence, due to the b — >• sj bound on the charged Higgs mass, this 
scenario is ruled out in the Type-II 2HDM. 



In the left panel of Fig. 11 we present the total i? 77 , with h, H and A summed, as a function of tan/3 for the 
constrained (red/black) and unconstrained (green/gray) scenarios in the Type-I 2HDM. This scenario does not differ 
much from the previous one where again only Type-I was allowed. In particular, when all constraints are imposed, 



the maximum value of R 11 has to lie very close to 1. In the right panel of Fig. 11 we present the allowed region in 
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Model I, h, H and A degenerate. Green - unconstrained; red - constrained Model I, h, H and A degenerate. Green - unconstrained; red - constrained 




1 2 3 4 5 0.5 1 1.5 2 2.5 

tan|3 r vbf 

YY 

FIG. 11: Left panel: Total R 11 (ft, H and A summed) as a function of tan/3 for the constrained (red/black) and unconstrained 
(green/gray) scenarios. Right panel: Allowed region in the R^ F — R 11 plane for the constrained (red/black) and unconstrained 
(green/gray) scenarios. 



the R^ ¥ -R 11 plane for the constrained (red/black) and unconstrained (green/gray) scenarios. If we now compare 
this scenario with the three previous ones we conclude that the differences are not that striking, especially in the 
constrained case: RY® F can be larger than the SM value reaching 1.7 while R 11 has to be very close to the value 
predicted by the SM. However, a large value of R^ F = 1.7 and i? 77 = 1 is only allowed in the present scenario 
and in the Type-I 2HDM with m# c± m^. In fact, the difference between these two scenarios is primarily due to 
the theoretical and experimental constraints: when tuh — rrih a nearly-fermiophobic scenario is allowed with values 
of R^ F reaching 3. In contrast, when all the neutral Higgs masses are degenerate the fermiophobic limit cannot be 
attained, and i?X? F reaches a maximum value of 1.7. 



8. CONCLUSIONS 

In this paper, we have addressed the possibility that an enhanced 77 signal in Higgs production at the LHC 
(under the assumption of a SM-like ZZ — >• 4 leptons signal) can be explained in the context of the 2HDM by a 
pair of nearly mass-degenerate neutral Higgs bosons with a mass around 125 GeV. To analyze this scenario, we have 
examined the softly broken Z2-symmetric and CP-conserving 2HDM with either Type-I or Type-II Higgs-fermion 
Yukawa couplings. We have scanned the resulting 2HDM parameter spaces, subject to the constraints of all known 
experimental observables (excluding the recently observed anomaly of B D^r~V T by the BaBar collaboration, 
which has not yet been confirmed by the Belle collaboration). 

We find that, in the degenerate mass scenario, it is only possible to produce an enhanced 77 signal in a region 
of parameter space where tan/3 is near 1. This result immediately implies that it is not possible to realize such a 
scenario in the MSSM, since the additional constraints imposed by supersymmetry combined with the LEP search 
for the Higgs bosons of the MSSM rule out values of tan/3 below about 2 [42] . 

In particular, we have demonstrated that an enhanced 77 signal could be due to a nearly mass-degenerate ft, A pair 
or ft, H pair. (In the case where H and A are nearly mass-degenerate, it is much more difficult to achieve a significant 
enhancement in the 77 signal subject to all the experimental constraints.) In the parameter region corresponding 
to an enhanced 77 signal, we generically expect an enhanced inclusive r + r~ signal, due to the contribution of the 
second mass-degenerate Higgs state in addition to the r + r~ signal produced by the SM-like Higgs state. 
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The absence of an enhanced r + r~ signal would be strong evidence against the scenario proposed in this paper. 
However, if both a 77 and r + r~ enhancement are confirmed, then one can probe the nature of the Higgs-fermion 
Yukawa couplings by distinguishing 77 signal events that arise from vector boson fusion (VBF). For example, we find 
that an enhanced VBF 77 rate is possible in Type-I models but is not possible in Type-II models. 

Of course, ultimately the verification of the mass-degenerate scenario would require the observation of two separate 
scalar states. As we do not expect these states to be exactly mass-degenerate, it is possible that the inherent 
experimental mass resolutions of the ATLAS and CMS experiments could eventually be sensitive to the presence of 
two (or more) approximately mass-degenerate scalar states. In particular, in the mass-degenerate scenarios studied 
in this paper, the 77 signal is a consequence of the production and decay of both mass-degenerate Higgs bosons, 
whereas the ZZ* — »• 4 lepton signal derives primarily from the production and decay of a single Higgs boson with 
SM-like couplings to vector boson pairs. Thus, an observation of slightly different invariant masses in the 77 and 
ZZ* — )> 4 lepton channels would be strong evidence that a scenario similar to the ones examined in this paper is 
realized in nature. 
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